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Abstract 

In this paper, we study the stabihty property of Hamiltonian systems on the Wasser- 
stein space. Let -ff be a given Hamiltonian satisfying certain properties. We regularize 
H using the Moreau-Yosida approximation and denote it by H-,-. We show that so- 
lutions of the Hamiltonian system for converge to a solution of the Hamiltonian 
system for H as t converges to zero. We provide sufficient conditions on H to carry 
out this process. 
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stability 

1 Introduction 

Let M. be the set of Borel probability measures on M.^ with finite second moments equipped 
with the Wasserstein metric. We study a Hamiltonian type evolution problem in A4. of 
following form : 



'1.1] 



where the given function H : M. ^ (— oo, oo] is referred to as a Hamiltonian. Here J : — )■ 
MP is a matrix satisfying Jvl-V for all v G M^. When D = 2d then we can simply set J to 
be the {2d) x (2c?) canonical symplectic matrix. Here, (9_if(/i) denotes the sub differential of 
H aX fi ^ Ai and T^Ai is the tangent space at /i in which will be defined below. There 
are various reasons for the terminology is Hamiltonian type. For example, fll.lj) is, roughly 
speaking, a limit of finite dimensional Hamiltonian ODE [8]. Geometric justification was 
also made in [9]. 
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The first systematic study addressing evolution problems on Ai of the Hamiltonian type 
was made by Ambrosio and Gangbo pQ. They studied the Hamiltonian system for locally 
subdifferentiable Hamiltonians and proved the existence of a solution. The theory in 
covers a large class of systems which have recently generated a lot of interest, including 
the Vlasov-Poisson in one space dimension the Vlasov-Monge-Ampere [5] [8] and the 

semigeostophic systems [3] [6] [7] [8] are casted into the Hamiltonian type formalism. 

We are interested in the stability property of Hamiltonian systems in the following sense. 
Let if be a given Hamiltonian. We ask ourselves whether there is any regular approximation 
Ht- of H such that solutions of (11. ip for H^. exist and converge to a solution of the system 
(11. ip for H as the approximation parameter r goes to zero. 

Since the Wasserstein space is an infinite dimensional metric space, the existence of such 
an approximation is not a simple question. In this paper, we show that the Moreau-Yosida 
approximation is the one we are looking for. Let if be a Hamiltonian satisfying assumptions 
(HI) and (H2) whose statements will be given later. We first regularize the Hamiltonian H 
to obtain Hr defined by 

HM = inf {l/2TW{fi, uf + H{u)]. 

The new functional Hr is 1/r— concave even if H is not. Next, we apply the algorithm 
developed in [T] to solve 

where d^Hr^jJil) is the superdifferential of at /i[ in the sense of [2]. Finally, we show, 
for any sequence r„ converging to zero, /i"^" (up to subsequence) converges to /i which is a 
solution of (II. ip . 

Our assumptions on the Hamiltonian H allow H to be no locally subdifferentiable. Hence, 
our stability result allow us to construct solutions to the system (II. ip for Hamiltonians which 
are not everywhere subdifferentiable around the initial measure. This is not the case in [1]. 
At the end of this paper, we will discuss more about how the Moreau-Yosida approximation 
scheme is useful in the study of non locally subdifferentiable Hamiltonians. 

We briefly summarize the contents of each section. Section 2 is a preliminary on the 
Wasserstein space M.. In section 3, we give an introduction to the Moreau-Yosida approxi- 
mation of functionals defined on Ai and investigate some properties of it. The main feature 
in this section is Lemma 13.41 which is the key ingredient to prove Theorem 14.51 In section 
4, we prove our main stability result Theorem 14.51 under assumptions (HI) and (H2) on 
the Hamiltonian H. We show Hamiltonians considered in [T] satisfy (HI) and (H2), and so 
corresponding Hamiltonian systems are stable w.r.t Moreau-Yosida approximation. Let us 
close this introduction by fixing notations and terminologies. 
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1.1 Notation and Terminology 

- V(R^) = is a Borel probability measure on M^} 

- Let Ai be the subspace of V{M.^) with bounded second moment, i.e. 

M := {fi eV{R^) : I2>0, dfi = l, \x\'^dfi<oo}. 

- Let fi G V{M.^) and let / : — )■ M'^ be a Borel map. Then z/ := is a Borel measure 
on M'^ characterized by I'lB] = fi[f~^{B)] for all Borel sets B C M'^. In this case, we say / 
pushes /i forward to v. 

- C^(M^) is the collection of all infinitely differentiable functions with compact support. 

- We denote Cb(M^) the collection of all continuous and bounded functions. 

- Let fin, fJ' £ P(M^), we define fin converges narrowly to fi if 

/ f{x)dfin{x) — > I f{x)dfi{x) as n — )■ oo, 

for any / G C{,(]R^), i.e. fin weak* converges to fi. 

- Id : — )■ is the identity map,z.e. Id{x) = x for all x G M"^. 

- 7r*,7r''-^ : R""^ — )■ R^,R^ x are the standard projections, i.e. 

'k'{xi,X2,- ■ ■ ,Xn) = Xi and 7r*'^(xi, X2, ■ ■ ■ , x^) = (xi, Xj). 
-Let fie V(R^) and let / : R^ ^ R^ We denote the norm of / by ||/||^, z.e. 

ll/IIJ:=ll/llW)= / \m\'dfi{x). 

- Let G P(R^), we denote the support of fi by supp{fi). 

- Let r > and a; G R'^ then Bx{r) denotes the open ball in R"^ of center x and radius r. 

- Let x,y e R^, we denote the inner product of x and y hj < x,y >. 



2 Wasserstein space 

Recall that Ai is the subspace of P(R^) with bounded second moment. In this section, we 
show that Ai has a metric structure and we introduce a differentiable structure in Ai. We 
refer to [2] and [13] for further details. 

2.1 Wasserstein distance 
Definition 2.1. Let /i, z/ G A^. Consider 

W2{fi,iy) := ( inf /" \x-y\^d-f{x,y)y^\ (2.1) 
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Here, r(/i, u) denotes the set of Borel measures 7 on IR-^ x which have /i and u as 
marginals, i.e. satisfying vr^(7) = /i and ii'^i'y) = v- 

Equation fl2.ip defines a metric on Ai which is called the Wasserstein distance. It is 
known that the infimum in the right hand side of equation (12. ip is always achieved. We will 
denote by To{fJ^, v) the set of 7 which achieve the minimum in (12. ip . 

Definition 2.2. Let /i, z/ G and 7 G ro(/i, v). The barycentric projection 7^ : — > M"^ 
of 7 with respect to the first marginal /x is characterized by 

^{x)%{x)dfi{x)= [ ij{x)yd^{x,y) y^j e (2.2) 



Similarly, the barycentric projection 7^ : M — )> M of 7 with respect to the second marginal 
z/ is defined by 

ij{yH{y)du{y)= f ^{y)xd-f{x,y) ^^j e C^. (2.3) 



2.2 Differential structure on Ai 

Definition 2.3. Given /i G A^, let T^A^ be the tangent space of at /i defined as the 
closure of VC^ in L'^{fi), i.e. 



T^M := {Vif : ip G C,°°(M^)} 

For any /i G Al, there is an orthogonal decomposition 

L2(/i) = T^M © [T^A^]^, (2.4) 

where [Tf,M]^ := {w G ^^(/i) : V ■ (w/i) = 0}. We will denote by tt^ : L'^ifi) T^M the 
corresponding orthogonal projection. 

As shown in [2], the tangent space enjoys many useful properties in analytic and geometric 
point of view. Here, we recall one of them which is related to absolutely continuous curves 
in Al. Let us first give the definition of absolutely continuous curves in metric spaces. 

Definition 2.4. Let (S, dist) be a metric space. A curve t G (a, 6) t— i- cr^ G S is 2-absolutely 
continuous if there exists /3 G L'^{a,b) such that 

dist(ai,a,) < [ f3{T)dT, (2.5) 

J s 

for all a < s < t < 6. We then write a G AC2{a,b;E>). For such curves the limit |cr'|(t) : = 
lim^^t dist((Tt, (Ts)/|t — s\ exists for £^-almost every t G {a,b). We call this limit the metric 
derivative of a at t. It satisfies \a'\ < (5 C^-a\m.osi everywhere. 

We now recall Theorem 8.3.1 in 0, which says that the tangent space provides a canonical 
velocity field for the absolutely continuous curves in Al. 
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Proposition 2.5. If fJ^ ^ AC2{a,b; Ai) then there exists a Borel map v : {a,b) x — )• 
such that Vt G L'^ifJ't) for -almost every t G (a, h) and 

^ + V ■ {v,^^,) = 0. 

We callv a velocity for fi. If w is another velocity for fi then 7r^i(t't) = 7Tm{wt) for -almost 
every t G (a, b), where tt^^ is defined in Definition \2.3[ Moreover, one can choose v such that 
Vt G T^^jvl and \\vt\\^^ = for C^-almost every t G {a,b). In that case, for C^-almost 

every t G {a,b), Vt is uniquely determined. We refer to v as the velocity of minimal norm, 
since if w is any other velocity associated to fi then \\vt\\^t < H^tH^^ for -almost every 
t G (a, 6) and so dist(/Xt, yUg) < Wv^W^j^^dr < WwrW^^dr for all a < s < t < b. 

Following [T], we give a notion of a differential and a definition of convex functions on 

M. 

Definition 2.6. Let H : M. ^ (—00,00] be a proper function on A^, i.e. the effective 
domain of H defined by D{H) := {/i G : H[fi) < 00} is not empty. We say that 
^ G L'^ifi) belongs to the suhdifferential d-H{fi) if 

H{u)>H{fi)+ sup / {^{x),y-x)d-f{x,y) + o{W2{fi,iy)), 

as u ^ fi. We denote the domain of suhdifferential by D(d^H) := {/i : d-H{fi) ^ 0}. 
If G d^{—H)(fi) then we say that ^ belongs to the superdifferential d~^H{fi). 

Remark 2.7. If d^H{fi) nd'^H{fi) 7^ then we say that H is differentiable at fi. In this case, 
there is a unique vector in d-H{fi) r\d^H{fi) (iT^Ai and we define the gradient vector V^iJ 
by the unique vector. 

Definition 2.8. Let H : M. ^ {—00, 00] be proper and let A G M. We say that H is A- 
convex if for every /io,/ii G and every optimal transport plan 7 G ro(/io,/ii) we have 

H{fit) < (1 - t)H{fio) + tH{fir) - ^t{l - t)l^|(/io, /ii) Vt G [0, 1], (2.6) 
where fit = ((1 - t)7r^ + t7r^)#l- If -H is (-A) -convex then H is called A- concave. 

3 Moreau-Yosida approximation 

In this section, we introduce the Moreau-Yosida approximation of functionals on Ai. 

Definition 3.1. Let H : M. ^ (—00,00] be a proper and coercive functional. For r > 0, 

the Moreau- Yosida approximation Hj^ of H is defined as 

HM= inl{^W^ifi,u) + Hiu)], (3.1) 
ueM I Zt ) 
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Here, H is coercive means that there exist r^, > and /i* G so that i/^-, (/^*) > —oc. We 
also set 

JM ■■= {f-^r : HM = + (3.2) 

Lemma 3.2. Let H : M. ^ (—00,00] he a proper and coercive functional, and Hr he the 
Moreau-Yosida approximation of H. Then Hr is ^ - concave. 

Proof. Let u E Aihe fixed, then it is well known that fi t— |W^2^(/i, z/) is a 1-concave function 
on M. This implies 

veM I Zt J 
is -—concave since it is an infimum of -—concave functionals. □ 

T T 

We now introduce two Lemmas which give relations between the sub differential of H 
and the super differential of Ht-. They play the key role in the convergence of Hamiltonian 
systems. 

Lemma 3.3. Let H : Ai ^ (—00,00] he a proper functional and Hr he the Moreau-Yosida 
approximation of H . For G given, if Vq ^ ^t[/^o] then Hr is superdifferentiahle at fio 
and H is suhdifferentiahle at Uo, i.e. /Xq G D[d^Hr) and G D[d-H). Furthermore, for any 
7 G ro(/io, 1^0), we have 

'^ed+Hrifio)nT^^M , ^ ed^H{iyo)nn^M (3.3) 

r r 
where 7^°(7^°y' is the harycentric projection of'~f with respect to the first (respectively, second) 
marginal as in Definition \2.^ 

Proof. From the definition of Ug G Jr[fio], we have 

Hrifi) - Hriflo) < ^Wiifi, uo) - ^Wiifio, ^o) G , (3.4) 
Zt Zt 

For a fixed we choose 77 G To{fio,f^)- Let 77 = Vxdfio{x) and 7 = lxdfio{x) be the 

disintegration of 77 and 7 w.r.t Ho- Define ui G V(M.^^) to be such that the disintegration of 
Ui w.r.t fio is 

/ VxX 'Jxdfioix). (3.5) 

We combine (13. 4p and (13. 5p to get 

Hril^) - Hr{llo) < T^W^K/^, ^o) - T^W^K/^o, ^o) 

Zt It 

1 /" \y — \x — 

< — / dw\[x,y,z) 

T Jj^sD 2 2 

= - / (x - 2;, y - x) H dui{x, y, z) 

T Jm?o z 

X — (x) 1 „ 

; y - x)d^{x, y) + -iy|(/Xo, /i) (3.6) 

r zr 
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which gives 



Id - 7,? 



T 

Furthermore, it is well known that Id — 7^° G T^^A^ (Proposition 4.3 of pj). This concludes 
the first inclusion of f l3.3p . 

To prove the second, we again exploit Uo G Jr[Aio] to get 

It It 
For a fixed let f] G Voiyo-, v) and define U2 G V{^^^ to be such that whose disintegration 
w.r.t Vo is 

73; X fj^duoix), 

where, fj = f^^^ fjxduo^x) and 7 = /jg^ 'jxdi^oix) are disintegrations of ry and 7 w.r.t Ug. 
Computations as in (13. 6p give 

Hr{iy) - Hr{Uo) > [ y - x)dfi{x, y) + 7^1^|(z/o, i^)- 

We conclude the second inclusion in (13. 3p using the same argument as above. □ 

Lemma 3.4. Let H : M. ^ {—00, 00] he a proper functional and let H^- be the Moreau- Yosida 
approximation of H . Given a sequence of measures /i„ and z/„ be such that z/„ G Jmlf^n]- 
Furthermore, suppose there is a constant C satisfying 

W2{fin,iyn) <CTn, (3.7) 

for all n. If fin converges narrowly to fi as Tn ^ 0, then z/„ also converges narrowly to /i. 
Furthermore, we have 

n c-o{V-^^^fir„ : n > m}) = n coH^^^^^^-^u^^ : n > m}). (3.8) 

m=l ^" m=l ^" 

i/ere cd denotes the closed convex hull with respect to weak* -topology. 

Proof. Narrow convergence of z/„ to yU is trivial from the assumption (13. 7p and the narrow 
convergence of /i„ to as r„ — > 0. 

To prove (EIH]), let us fix G C^{M.^- M.^). Then we have 

{^{y).^-^^^^{y))du^M= I {Hy),^)dirAx,y) 

X — y 

{i/j{x) + V^(^^,y) ■ {y - x), )d-fr^{x,y) 

(^(x), (x))rf/i„(x) 

Tn 

+ / W(e.J-(2/-a:),^^)ci7r„(x,t/) (3.9) 
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where 7^-^ G Toifim^^Tn) and ^x,y is a point on the hne segment connecting x and y. Since 



oomD. Tn>D\ 



we have 



x-y, 



We combine equations fl3.9p and f l3.10p to get 



< 



liv^ll 



< IIVV^IIooCV 



(3.10) 



hm / (?/))c^^.„(2/) = hm / (^(x), 



which concludes fl3.8p . 



□ 



4 Convergence of Hamiltonian systems w.r.t Moreau- 
Yosida approximation 

Now we are ready to state our main result on the stability of Hamiltonian systems. More 
specifically, solutions of the approximated Hamiltonians converge to a solution of the original 
Hamiltonian system. Let us first be clear about the meaning of solution. 

Definition 4.1. Let H : Ai ^ (—00,00] be a proper and lower semicontinuous function. 
We say that an 2— absolutely continuous curve fit '■ [0,T] — )■ D{H) is a solution of the 
Hamiltonian system starting from /2 G A^, if there exist a vector field Vt G L^ifit) with 
||ft||^t G L^(0,T), such that 



|yUt + V ■ (Mm) = 0, fio = fi tE{0, T) 
Vt G d_H{fit) n T^^M a.e t G (0,T). 



(4.1) 



Equation fl4.ll) should be understood in the sense of a distribution: For any 77 G C^(0,T) 
and C G CI 



;~(M^ 



), we have 

T 







r]\t)C{x)+r]{t){VC{x) : Jt;t(x))rf/Xf(x)rft = 0. 



To ensure the stability of Hamiltonian systems, we require two assumptions on the Hamil- 
tonian. 

(HI) There exist constants Co G [0, 00), -Ro G (0, 00] such that if W2{fi,Ji) < Ro then 
fi G D{H) and , for each fi, there exists a unique v G t/r[/w] satisfying 

• fi ^ u E Jrlfi] is continuous w.r.t the topology induced by the Wasserstein distance 
and 

< C„. (4.2) 
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• There exists a constant k > such that 

If supp{fi) C i?o(^) then supp{u) C BQ^kr), (4.3) 
for all r > and fi. Recall, -Bo(^) is the ball around the origin with radius r in M^. 
(H2) If fin £ D{d^H) and converges narrowly to /i, then /i G D{d^H) and we have 

oo 

Pi co{{Wnlin ■■ Wn G d_H{fin) H T^,M, u > m}) C {w fi : w G d_H{ii) n T^A^}. (4.4) 

m=l 

Remark 4.2. 1. Notice that our Hamiltonian if does not need to be subdifferentiable 
everywhere in a neighborhood of Ji. We only assume that D[d-H) is closed in the 
weak* topology and (14.41) holds. 

2. Concerning (HI), suppose H satisfies the following convexity condition for some A G M: 
For all yU, z/q and ui in D{H), there exists a curve cr : [0, 1] — such that ao = h'o,ai = 



ui and 



1 + 



nir, /i; at) < (1 - t)'H(r, /i; z/q) + t?^(r, /i; z/i) t(l - t)W^{iyo, i^i), (4.5) 

2r 

for all t G [0, 1] and < r < ^. Here, H(r, /i; z/) := ^Wi{iy, /i) + ii(z/). 

Then, Theorem 4.1.2 in [2] says that if /i G D{H) and Ar > — 1 then there exists a 
unique Hr G JrM and the map n G D{H) fir ^ JAfA is continuous. 

4.1 Existence of solutions of the regularized Hamiltonian systems 

Lemma 4.3. Let H : Ai (— oo, oo] be proper and lower semicontinuous, and satisfy (HI ). 
Let fin be a sequence satisfying W2{fin,Jj) < Ro and z/„ G Jr[fin]- If fJ^n converges to fi in the 
Wasserstein topology, then z/„ also converges to u E JrifA in the same topology. Furthermore, 
we have 

ICo := fl co({ -^fin : n > m}) C {— ^y^^}, (4.6) 

m=l 

where 7„ G To{fin, ^n) and 7 G To{fi, v). 

Proof. By (HI), there exists a z/ G Jt[iA such that W-zi^n^ ^) — ^ 0. 

Next, to prove (14.61) . let us assume u G /Co- Then, there exists a sequence {Am}m=i such 
that 

A„^ = 5ZA^ 5^Ar = l, 0<Ar<l, m</^GN 
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and Km weak* converges to u. For any F G Cc(M''^; M-^), we have 



[ F ■ ciu = lim V / {F{x) 



T 



i=m 
I 



limVAr/ (F(x),^)rf7.(a:,|/), (4.7) 



where 7, G To{fii, i^i)- Since W^2(/^n., W2{i^n, z^) — as — 00, there exists 7 G To{fi, u) so 
that 

hm W^2(7i,7) = 0. (4.8) 

We combine (HTj) and (CSj) . to get 

F-dM= I {F{x),^)d^{x,y), (4.9) 



which proves (14. 6p . □ 

Now we generate a solution of the Hamiltonian system for Hr- The proof of the following 
theorem is based on Theorem 7.4 in [1]. 

Theorem 4.4. Let H : M. (— C)0, 00] he a proper and lower semicontinuous functional 
satisfying the assumption (HI). Let Co and Ro be constants in (HI), and set T = 

If pL E M. has hounded support, then for each r > 0, there exists a solution of the following 
Hamiltonian system starting from fl 

+ V ■ (K/i[) = 0, /i5 = /i tG(0,T) 

vl = G d+Hr{fil) n T^jM, a.e t G (0, T), 

where vl G Jr[^l]. Furthermore, t ^ fil is Lipschitz continuous w.r.t the Wasserstein dis- 
tance with Lipschitz constant Co which is independent of r. 

Proof. Step 1. Construction of a discrete solution 

For given r > 0, we fix an integer and divide [0, T] in equal subintervals of length 
h = T/N. We build discrete solutions satisfying: 

(a) The Lipschitz constant of t 1— )■ /if!^ G is less than Co- 

(b) For all t G [0,T], we have supp{^f^) C -Bo(e r) if supplfi) C BQ{r). Here, k > is 
same as in (HI). 

(c) /if!^ satisfies the discrete Hamiltonian equation 

l/^f. + V ■«/.?;) = 0, tG(0,T), (4.11) 
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with 



Id - 7 AT 

u;f^ = J ^ for t = 0,h,2h,--- ,Nh, (4.12) 

T 

where u^^ G Jr[^fr] and 7 G ro(/if^, u^^). 

Since and r are fixed, we use the notation fit := fif^^ for convenience, 
(i) We build the solution in [0, h]. 

Let us call /xq := and choose uq G JtI/Uo] by (HI). We fix 7 G To{fiQ, 1^0) and set 



Id - 7^^° 

Wo := i 

We define 



Wo:=I ^. (4.13) 



/ . , , , X (/c; + two)#(wo/io) , ^ rn Li 
/it := [Id + two)#/io, := , t G [0, h\. 

We claim that Wt is a velocity field for fit, that is, 

+ V ■ (M;t/iJ = 0, (4.14) 

dt 

holds in the distribution sense in (0, h). Indeed, for any G C^(M^), we have 



— / (f)dfit = -77 / 0(-^'^ + two)dfio 
at j^D at J^D 

{V(j){x + two{x)),Wo{x))dfio{x) = I {V(f),wt)dfit. (4.15) 
Notice that Lemma 7.1 in [T] gives 

\wt\'^dfit < I \wo\'^dfio, (4.16) 



for all t G [0, h]. Jensen's inequality with (HI) gives 



\wo\'dfio= ^dfio<-W,'{fio,iyo)<C',. (4.17) 

We exploit Proposition 12.51 with fl4.16p and fl4.17p . to conclude that t 1— /ij is Lipschitz 

continuous with a Lipschitz constant Co- 
Next we show the bound on the support. Since supp{fio) C Bq^t), we have suppi^uo) C 

Bo^kr) by (HI). Hence, if 2; G supp{fit) then 

I ~ 7un(^) I / , r + kr. 
\z\< sup |x + a ^^^\<{r + h— — ), 

x£supp{po) "7" 
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hence we have supp{fit) C -Bo((l + ^^^)^)- 

(ii) We continue this process in [h,2h]. 

Since we have W2{fio, fih) < hCo < Ro, we can choose G JrlfJ'h] and set 

Id - 7,';" 

Wh:=I ^. (4.18) 

r 

We introduce the following extension for t G {h, 2h], 

[id + (t - h)wh)Jwhfih) 
fit = {Id+ (t - h)wh)JfJ^h), wt = . 

As in (i), we can check t i— )■ /i^ is Lipschitz continuous with a Lipschitz constant Co in 
[h, 2h]. Furthermore, the equation fl4.15p holds and we have supp{fit) C Bq{{1 + ^^hYr) for 
allte[h,2h]. 

(iii) We iterate the above process until we get a Lipschitz curve t y-^ fit E M. with Lipschitz 
constant Co- The curve satisfies fl4.15p for a.e t G (0,T) and hence 

d ^ , s 
—fit + V ■ [wtfitj = 0, 

holds in the distribution sense. Furthermore, for all t G [0, T] 

suppifxt) C 5o((l + — fr) C Bo{e^—r). 

Recalling that /if!^ := fit completes the proof. 

Step 2. Let increase to oo 

From (a), t h-> fi^^ are equi-bounded in Ai and equi-Lipschitz with Lipschitz constant 
Co- Since fif^^ have uniformly bounded supports, we may assume(up to a subsequence) that 
fif^^ converges in the Wasserstein topology as — > oo. That is, there exist /i[ such that 
W2{fif^^, fil) — for any t G [0,T]. Moreover, t H- /i[ is Lipschitz continuous with Lipschitz 
constant Co- As shown in [1], /xj" solves + V ■ [wlfi].) = with the following property 

oo 

wlfil G fl co{<Xr --NyM} a.e t G (0, T). 

M=l 

Since 

w^^fil = {Id+{t - [Art]/Ar)i/;0^,]/^_J^(«;5^,]/^_,^f)^,]/^_J 

= {Id+{t - m/N)^ ^^^^^^tlN,r)■ 



We also obtain 



id-^T""-^ 

e n '-^^^l^mN,r : A^ > M}. (4.19) 

M=l 
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Lemma [4.31 together with (14.191) gives 

Id - fJ 

T 

where G Jrif^l] for a.e t G (0,T). This with Lemma [3.31 concludes the proof. □ 



4.2 Stability of Hamiltonian flows 

Theorem 4.5. Let H : M. ^ (—00,00] he a proper and lower semi- continuous functional 
satisfying (HI ) and (H2). We assume that Ji ^ M. has a hounded support. For each r G (0, 1), 
let jX^ he the solution of the system (^T7^ in Theorem \4-4\ Then, {/i'^}T>o (up to a sequence) 
converges to a solution of the Hamiltonian system 

ifit + V- (M^t) =0, fio = fi t G (0, T) 

vt G d_H{fit) n T^,M a.e t e (0,T), ^ ^ 

as T converges to 0. 

Proof. Since t ^ nl are equi-bounded in M. and equi-Lipschitz continuous, we may assume 
that, for any t G [0,T], converges narrowly to Ht as r„ — >■ 0, for some subsequence r„. 
By the same reasoning as step 2 in the proof of Theorem 14.41 /xj solves 



+ V ■ (Ji;t/it) = 0, tG(0,T) 



(4.21) 



with 



for a.e t G (0,T). Here 



00 



Vt^Jit e Pi co{t;["/i[" : n > M}, (4.22) 



M=l 



Id — ■y'^%„ 

vl- = ^ G 9+i7(/i-) n T^^.M, vl- G J.„ [/i-]. (4.23) 

From Lemma WM together with (14.221) and (I4.23p . we know vl"^ — /i^ narrowly as r„ — )■ and 

00 

vt[it e fl co{e["i^r : n > M}, (4.24) 

M=l 

where z/[" G Jr,J/i[1 and = — G d^HiuJ") n T^^nAl. 
By (H2) and flCTj) . we get 



i;* G 9_i7(/it) n T^,7W, a.e t G (0, T), 
which concludes the proof. 



□ 
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4.3 Example 

Let Ho & M. have a bounded support and a > 0, we define 

H{fi) := + J V{x)dfi{x) + JJ W{x,y)d^I{x)d^I{y), (4.25) 

where V : — M is Ay— convex for some Ay G M, and VI^ : x — )■ M is convex and 
even. Assume also that both are differentiable and have at most quadratic growth at infinity. 
Then, the function H : M. ^ (— C)0, oo] as in fl4.25p satisfies (HI) and (H2). 

Proof. First, we notice that H defined as ()4.25p is (Ay — a)— convex and locally Lipschitz(pP), 
i.e. there exist i?^, > such that 

i/(/ii) - i7(/i2) < C^1^2(/ii, /i2), (4.26) 

for all Hi with W2{fi, fn) < i = 1,2. 

Secondly, H satisfies the convexity condition (14. 5 p with A = Ay — a (Chapter 9 in [2]). 
From Remark 14.21 it follows that for all sufficiently small r > and V/i G Ai, there exists a 
unique fir G Jrifj]- Furthermore, for fixed ji ^ jir ^ JAlA is continuous. 

We now show there is Q < Ro < Rfi such that for all sufficiently small r > 0, 

if W2{H,fL)<Ro then W2{ft, ^^r) < Rfi. (4.27) 
Once we have f l4.27p . (14.261) with fir G JAfA gives 

for all yU with W2{H:P') < -^o- That is, 

-Vr2(/i,/i.) <2C^, (4.28) 
r 

which implies that (jO]) in (HI) holds with Co = 2Cfi. 

Now, let us prove fl4.27p . We define Ro := Rp,/2. If H^2(At, /i) < Ro then we have 

^W-2'(/^,/^r) + i^(/X.)<i/(/x) 

<H{fi) + CpW2m,fi) 

< H{fi) + CpRo := C. (4.29) 

We need to estimate H{nr) in (14.291) . Since V and have at most quadratic growth at 
infinity, there are constants ci and C2 such that 

|^(a;)| < ci|x|2 + C2 and \W{x,y)\ < Ci{\x\^ + \y\^) + C2, (4.30) 
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for all x,y G 'MP . This implies 



a. 



> -^W^2^(/^r, /io) - 3Ci j \x\^dfir{x) - 2C2 



= -^Wi{fir,f^o) - 3ciiy|(/i.,5o) - 2C2. (4.31) 
We combine f l4.29p and f!4.3ip . and get 

;^l^22(/i,/i.) < ^W^ifir,fio) + Sc,W^ifir,6o) + 2c2 + C. (4.32) 

Now let r —J- in (14.321) . Since /io, 6o are fixed and lV2(/i5 fJ') ^ Ro = -R/t/2, 

W^2(/^5 /^t) — ^ uniformly w.r.t /i as r — 0, (4.33) 
which implies f l4.27p . 

To finish proving (HI), it remains to prove (14. 3p . Let H : M. ^ (— C)0, oo] be defined by 

For the Hamiltonian H, it was shown that (14.31) holds for some A; > [refer |10]). It is easy 
to see (14.31) holds with same /c > for the Hamiltonian H as in (I4.25p . 

Hence, for the Hamiltonian H defined by (I4.25p . the assumption (HI) holds with Co = 
2Cfi and Ro = Rp./'i. It was shown in [1] that the assumption (H2) also holds. 

□ 

Comments: Suppose we want to solve the finite dimensional Hamiltonian system which 
consists of a single particle 

- x"(t) = -VV{x{t)) 
x'{^) = v,x{^) = x, ^ ^ 

where V : MP — )■ M is given and -y, x G M^ . If V is not everywhere differentiable then we 
may try a regularization scheme as follows. We first solve the approximate system 

x:'(t) = -vK(x,(t)) 

x^(0) = ti,a;,(0) = x, ^ ' 

where K is any regular approximation of V. For example, we can define := pe *V a,s 
the standard mollification of V . Next we check if the solution x^{t) of (14.351) converges to 
a solution x{t) of (I4.34p as e goes to zero. Of course, we need certain properties on V to 
ensure this stability property hold. For example, if V is convex then the limiting solution 
x{t) satisfies the differential inclusion x"{t) G —d-V{x{t)) instead of the first equation in 

diii. 
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Let us address the Hamiltonian system (11. ip in the Wasserstein space. As we saw in the 
previous sections, under certain conditions on the Hamiltonian if, the Hamihonian system 
is stable with respect to the Moreau-Yosida approximation. Therefore, we may apply the 
Moreau-Yosida approximation scheme to study non locally subdifferentiable Hamiltonians. 

Let p, E Ai and i3 be a neighborhood of ft in the Wasserstein space. Suppose our 
Hamiltonian H is subdifferentiable only in a proper subset V G B, and p, E T). We want 
to solve the system (11. ip with the initial measure Jx. To do this, we need an algorithm to 
construct solutions which stay in the subset V. 

In Theorem 14. 4[ we construct approximate solutions for H as well as solutions /i^ for 
Hr- Notice that we have Ur E V and we need only the assumption (HI) on H in Theorem 
14.41 Note, the assumption (HI) has nothing to do with the subdifferentiability of H, which 
means that the construction of approximate solutions Ur relies entirely on the Moreau-Yosida 
approximation method. Next, in Theorem 14. 5^ we add the assumption (H2) on H which then 
implies the convergence of Ur to a solution of the system (II. ip . The assumption (H2) does 
not require our Hamiltonian H to be subdifferentiable everywhere in the neighborhood B of 
p,. Instead, it requires V is closed in the weak* topology and (14. 4p hold. 

Hence, as a direct result of our stability result, the Moreau-Yosida approximation scheme 
provides an algorithm to construct a solution of the system (II. ip for Hamiltonians which are 
subdifferentiable only in a proper subset D of a neighborhood of p. 

Acknowledgements : I would like to acknowledge my gratitude to Jacob Bedrossian 
and Ryan Hynd for valuable suggestions. 
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